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ABSTRACT 

We present a reanalysis of the X-ray data for RX J0720. 4-3125 presented in Zane et 
al. (2002), using more data recently available from XMM-Newton and Chandra. This 
analysis also corrects the ROSAT data used in that paper to the TDB time system, 
incorporates the revised XMM-Newton barycentric correction available since then, 
and corrects the definition of the instantaneous period in the maximum likelihood 
periodogramme search. However, we are now unable to find a single coherent period 
that is consistent with all ROSAT, Chandra and XMM-Newton datasets. From an 
analysis of the separate datasets, we have derived limits on the period change of 
P = (1.4 ± 0.6) x 10~ 13 s/s at 99% confidence level. This is stronger than the value 
presented in Zane et al. (2002), but sufficiently similar that their scientific conclusions 
remain unchanged. We examine the implications in more detail, and find that RX 
J0720. 4-3125 can have been born as a magnetar provided that it has a young age of 
~ 10 4 yr. A more conservative interpretation is that the field strength has remained 
relatively unchanged at just over 10 13 G, over the ~ 10 6 yr lifetime of the star. 
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1 INTRODUCTION 

RXJ0720.4-3125 is a nearby, isolated neutron star (INS) 
originally discovered with ROSAT during a systematic sur- 
vey of the Galactic plane by Haberl et al. (1997). A clear 
modulation of the X-ray intensity is detected at a period 
of 8.391 s, which identifies the spin velocity of the neutron 
star. 

RX J0720. 4-3125 is, so far, the isolated neutron star 
with best studied timing properties. Its proximity and rel- 
atively high brightness made possible not only to measure 
the spin period, but also to study the period changes over 
a long-term. This is crucial, since the positive detection of 
spin-up or spin-down phases can shed light on the mecha- 
nisms that regulate the neutron star interaction with its sur- 
roundings. A large positive spin-down can indicate magneto- 
dipolar breaking, leading in turn to an estimate of magnetic 
field strength and star spin down age. In the attempt to 
measure a secular period change, Zane et al. (2002, here- 
after Z02) recently undertook a comprehensive timing analy- 
sis of ROSAT, BeppoSAX, Chandra and XMM-Newton data 
spanning a total period of ~ 7 yrs. Independently, a similar 
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analysis based on an additional Chandra dataset, but with- 
out the XMM-Newton data, was presented by Kaplan et 
al. (2002, hereafter K02). Although the two studies agreed 
in their main conclusion that the period change was less 
than ~ 3 x 10~ 13 s/s, it was evident that there were in- 
consistencies in the details of the analysis which have impli- 
cations for further studies. The arrival times computed by 
the two groups were different. K02 also raised theoretical 
concerns about the validity of the coherent analysis of the 
entire dataset carried out by Z02. That was the only method 
that, at that time, permitted the (positive) sign of the pe- 
riod change to be constrained. In fact, the phase-incoherent 
analysis presented by Kaplan et al. (2002) still leads to a 
large uncertainty in P, and does not permit the discrimina- 
tion between positive and negative spin-down. 

The above-mentioned inconsistencies have motivated us 
to undertake a complete revision of the solution. At the time 
Z02 was published, timing inconsistencies between XMM- 
Newton and radio data on other pulsars had come to light 
(Kuster et al. 2001). This was traced to an error in the 
XMM-Newton Science Analysis System (SAS) task barycen 
before version 13.1 [the spacecraft position vector was con- 
structed as (x, y, y) instead of (x, y, z)]. Although we argued 
in Z02 that any effect on the coherent analysis was small, 
we nevertheless performed a reanalysis of the XMM-Newton 
data to check, and also to investigate whether this corrected 
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the inconsistencies with the K02 analysis. We found indeed 
that the effect was negligible. The implication is therefore 
that the inconsistencies between the Z02 and K02 analyses 
remained. 

We have therefore undertaken a complete reanalysis of 
the dataset to investigate and eliminate these inconsisten- 
cies. This has brought to light that the times used for the 
ROSAT data in both the Z02 and K02 studies are in UTC 
rather than TDB as used for Chandra and XMM-Newton; 
the times of arrival (ToAs) in K02 are incorrectly calculated; 
and that the formula used to calculate the instantaneous 
phase was incorrect in Z02. Substantial new datasets are 
also available from both XMM-Newton and Chandra which 
permit a more comprehensive investigation into the period 
evolution. We report here the results of this reanalysis. 



2 NEW OBSERVATIONS AND REVISIONS 

We refer the reader to Z02 and K02 for the details of the 
previous observations, and comment here largely only on 
those aspects which have changed. The table of observations 
is given in Table Q 



2.1 XMM-Newton Data 

RX J0720.4-3125 has observed four times by XMM-Newton, 
in revolutions 78, 175, 533/534 and 622. This is a substantial 
increase over that available to Z02. 

In Z02, for revolution 78 (XOOa in Z02), special arrange- 
ments (tcs_fix) were made to ensure the timing accuracy of 
the data in the absence of a full observation data file (ODF). 
In the intervening period, the observation has been repro- 
cessed in the XMM-Newton pipeline, so that the standard 
observation data file (ODF) is now available. We used SAS 

5.2 odfingest to prepare the ODFs and were also assisted by 
U. Lammers at ESA who ran a more recent (not then pub- 
lic) version. Nothing unexpected with regard to timing could 
be found in the logs or in a direct inspection of the ODF 
time correlation and orbit data. In running the pipelines 
(SAS emchain and epchairi), we tracked down timing cor- 
rection extrapolation warnings in EPIC-pn to anomalous 
auxiliary data within the first 100 frames after the expo- 
sure start in some quadrants, causing the computation of 
the time tags to fail. This problem was fixed by U. Lammers 
using software which was later released in SAS 5.3.2. Also 
in EPIC-pn for the XOOb data, we corrected the 1 s errors in 
the frametime as described in Z02. For the EPIC-mos data 
we found some timing warnings referring to non-increasing 
times: these were flagged automatically and not included in 
the later selection procedure we used to create the extracted 
event lists. In other regards, the extraction was performed 
as in Z02. 

The rev. 533/534 observations of RX J0720.4-3125 were 
made for calibration purposes on 2002 Nov 6-9. Extending 
over 3 days, these data can provide tighter constraints on 
the period determinations than the earlier observations. The 
EPIC-pn data only were extracted using the procedures de- 
scribed in Z02, with the same extraction parameters. SAS 
v5.3.2 was used to process the events. Corrections were ap- 
plied for the 1 s frametime errors as described in Z02. Rev 



622 observations were also made for calibration purposes on 
2003 May 2. These were processed as for rev 533/534. 

For both old and new datasets the barycentric correc- 
tion was carried out using SAS barycen version 13.1 and 
13.2. 

2.2 ROSAT Data 

Despite that the ROSAT data had been processed in two dif- 
ferent ways in Z02 (as a double check, using both the ROSAT 
EXSAS software system and FTOOLs), an error was made 
in not converting the time system from UTC to TDB in or- 
der to relate these observations to the XMM-Newton and 
Chandra data, both of which use the TT/TDB system. The 
timings for the ROSAT data were therefore recomputed, us- 
ing the convert /utc_tdb task in EXSAS. The absolute timing 
of the 1998 ROSAT data is uncertain, which prevents their 
use in any coherent analysis (Z02). 

We have also extracted the data for RX J0720.4-3125 
from the ROSAT All Sky Survey. These data suffer from 
relatively low count rates because of vignetting away from 
the field centre, and also the obscuration from the PSPC de- 
tector support grids. On the other hand the duration of the 
observations is several days (2.07 d) days, so they potentially 
provide a good constraint on the period determination. Fi- 
nally we have also extracted early HRI observations taken in 
1994 May. The source is near the edge of the detector area, 
and the observation yields only a small number of counts. 
Nevertheless, as for the RASS data, the duration of the ob- 
servation provides some useful constraints. 

2.3 BeppoSAX Data 

In Z02, the BeppoSAX observations were not directly used 
in the period analysis, but were folded on the derived pe- 
riods from the remaining data to check the consistency of 
the period solution. As the FTOOLS Earth2Sun correction 
used in Z02 did not apply the BeppoSAX clock correction, 
the BeppoSAX data have been reprocessed using the SAX 
Data Analysis System and the correction recomputed using 
the baryconv task in SAX DAS (kindly performed by M. 
Feroci). It was confirmed by T. Oosterbroek (private com- 
munication) that the SAX timings after this operation are 
in the TDB system. 

2.4 Chandra Data 

Since the work of Z02, Chandra observations made in 2001 
December 4-6 have become publically available. These ob- 
servations were reported in K02. In this work, we used the 
Chandra data on RX J0720. 4-3125 collected in four obser- 
vations of 2001 December 4-6, with effective exposures of 
15.0, 10.6, 4.1 and 1.9 ks. The observations were conducted 
with the Advanced CCD Imaging Spectrometer (ACIS-S) 
operated in Continuous Clocking mode. This observational 
mode provides a 2.9 ms time resolution by means of sacri- 
ficing spatial resolution in one dimension. In total, we used 
66,822 counts extracted from segments of a 3'.'44 width in 
the one-dimentional images of the source, in the 0.1-1.0 keV 
range. The times of arrival were corrected for the dither and 
the Science Instrument Module motion as described in detail 
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Date 


Observatory 


Instrument 


Exposure 


Exposure 


Effective 


Label 








identification 


duration 


exposure 












(s) 


(b) 




1990 Oct 11 


ROSAT 


PSPC 




178840 




RASS 


1993 Spn 27 


ROSAT 


PSPC 


rr>300338n00 


11980 


3221 


R93 


1 994 Mav 1 1 


TIGS AT 


HRI 


rh201 733n00 


173237 


7402 


R94 


199fi Anr 25 


TIGS AT 


HRI 


rri300508n00 


7743 


3125 


R96a 


1 QQfl Mav 7 


JiftQA T 


HRI 


1 ll_LO<J 1UUI1UU 


7838 


3566 


R96b 


199fi Spn 27 


ROSA T 


HRI 


rh300508a01 


1498 


1409 


R96c 


199fi Nnv 3 


ROSA T 


HRI 


rh400884n00 


65698 


33569 


R96d 


1 997 Mar 1 fi 




LECS 


T.ECS 20979091 


99418 


17235 


S97 


1998 Apr 20 


ROSAT 


HRI 


rh400944n00 


460195 


3566 


R98 


2000 Feb 1 


Chandra 


HRC-S(LETG 1st order) 


348+349+745 


305528 


37635 


ChOO 


2000 May 13 


XMM-Newton 


MOS1 + thin filter 


0124100101-001 


61352 


61648 


XOOa 






MOS2 + thin filter 


0124100101-002 


61648 


61648 








PN + thin filter 


0124100101-003 


52305 


52305 




2000 Nov 21 


XMM-Newton 


MOS1 + medium filter 


0132520301-007 


17997 


17997 


XOOb 






MOS2 + medium filter 


0132520301-008 


17994 


17994 








PN + medium filter 


0132520301-003 


25651 


25651 




2001 Dec 4 


Chandra 


ACIS-S 


2771-2774 


171243 


31532 


ChOl 


2002 Nov 6 


XMM-Newton 


PN + thin filter 


0156960201-003, 0156960401-003 


208582 


58554 


X02 


2003 May 2 


XMM-Newton 


PN + thick filter 


0158360201-023 


72793 


72793 


X03 



Table 1. The ROSAT, BeppoSAX, Chandra and XMM-Newton observations of RXJ0720.4-3125 used in this paper. 



cp_new_chandra_2.ml.ps (2003 Apr 1 i at 13:15:39) cp_new_xmm.ml.ps (2002 Dec 20 at 10:27:50) cp_ncw_xmm_6.ml.ps (2003 Jul 20 at 13:15:05) 




8.3900 8.3905 8.3910 8.3915 8.3920 8.386 8.388 8.390 8.392 8.394 8.3900 8.3905 8.3910 8.3915 8.3920 



Figure 1. The MLPs for the ChOl data (left), the X02 data (center) and the X03 data (right). The 68 percent ltx confidence levels are 
set from Ax 2 = 1, and are not discernible on these plots. 



in Zavlin et al. (2000), and transformed to the Solar System 
Barycenter using the axBary tool of the CIAO package 1 . 

3 MODEL REVISIONS 

Z02 (equation A5) defined the instantaneous period to be 
P(t) = P + P(t - to) , (1) 



http : //asc .harvard. edu/ciao/ 



where Po is the period at t = to and P is the period change. 
This period was substituted in the model prediction used in 
the Maximum Likelihood Periodogramme (MLP) (equation 
A4 in Z02) 

/(*) = oo[l + Acos(2wt/P(t) + 6a)} , (2) 

where A is the normalised amplitude and 8o is the phase. 

This substitution is a rough approximation: as the co- 
sine argument in equation J3J is a phase corresponding to 
the orientation of the neutron star, the cumulative outcome 
of the past history of the period variations is not negligible. 
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w(t) 



The definition of angular frequency is: 
2tt d(/> 



P(t) dt ' 

thus, substituting P(t) from equation Q, the phase at time 
t is 

rt 



4>(t) - 0(t o ) = 2tt 



1 



P - Pto + Pt 



-dt. 



(3) 



In our prescription, P is constant, so this can be integrated 
analytically to 



0(t)-0(i o ) = ^ln^l 
P P) 



(4) 



Accordingly, the model prescription should be revised from 
equation to 



lit) = a 



1 + Acos(^ln^+0 O ) 
P Po 



(•>) 



where P(t) is as defined in equation (y and all of the timing 
and phase zero points are included in Oq. 

The calculation of the phase in the MLP via equa- 
tion © is slower than that previously done via equation 
and care is required in order to minimise the numerical errors 
resulting from the calculation of the logarithm of a number 
close to unity. In particular in the case where P = 0, the 
solution of equation (3) is <j>(£) — <j>(to) = 2n(t — io)/P) and 
then the intensity model is simply 

I(t) = ao[l + Acos(^+6o)]. (6) 
Po 

It could be argued that in the presence of constant 
torques, a constant !> assumption is more valid than a con- 
stant P. In this case the accumulated phase is <f>(t) — </>(to) = 
2nv(t) + 2-Kvt 2 /2. In practise we checked that, for P = 
5 x 10~ 13 s/s (Z02, K02), the accumulated phase difference 
over ten years between these assumptions is < 0.5 percent 
and therefore negligible. The two assumptions give consis- 
tently similar results, and this also implies that a small sec- 
ond derivative term v (of the kind v — 2P/P 3 — P/P 2 ) does 
not affect the solution. With the increase in the time base- 
line, however, it will become increasingly possible to distin- 
guish between the models. In order to facilitate comparison 
with Z02 and K02 we retain the constant P model here. 



step- wise fashion. First the best fitting periods and their un- 
certainties have been derived for the new Chandra (ChOl) 
and XMM-Newton (X02, X03) datasets (FigureQ, and also 
the S97 and the RASS datasets using a MLP as in Z02. 
These were added to those derived in Z02. The full list is 
given in Table The first stage is to use the XMM-Newton, 
ChOl, R93 and R96d datasets to determine unequivocally 
the appropriate alias peaks in the ChOO data and the R98 
data, since in these individual datasets there is no question 
of the period determination. We performed a weighted least 
squares fit using these data, with parameters Po and P and 
the extrema of the 99 percent confidence ellipsoids deter- 
mined. The result is shown in Figure [5] (top). The solid line 
is the best fit P = Po + Pt, while the dotted lines bound 
the range in P permitted from the 99 percent extrema. Also 
shown are the aliases of the period determinations for the 
R96 and ChOO data. It is clear that only the aliases cor- 
responding to P ~ 8.3911 s are consistent with the least 
squares fit. These indeed correspond to those selected in 
Z02 and K02, but this (at least in Z02) relied on the corre- 
spondence of these peaks in the R98 and ChOO MLPs. 

It is now possible to include the periods for these peaks 
in the least squares fit and repeat the process. The result 
is shown in Figure |2] (middle) . Also shown are the aliases 
of the combined R96a and R96b determinations: these are 
short runs but have a convenient two- week separation. Again 
from the 99 percent confidence range, the central alias at 
P ~ 8.3911 s is the only acceptable period. The 99 percent 
confidence interval changes slightly if we include the RASS 
period determination, but the selection of the aliases is un- 
affected. This in turn can be used to further constrain the 
fit, and to determine which alias should be selected from 
the combined R96c and R96d data, which have a 1 month 
separation. This is shown in Figure [5] (bottom), where the 
99 percent confidence limits identify the alias at 8.391095 s 
as the appropriate period. Again, the RASS data can be in- 
cluded without any change to the alias selection. Including 
the period determinated from this alias, the least confidence 
intervals in the Po , P plane are shown in Figure |H] 

The outcome of this analysis is that the 99 percent con- 
fidence interval for the P term lies in the range 0.70 x 10~ 13 
to 2.10 x 10 -13 s s _1 , with the 68 percent interval being 
1.05 x 10" 13 to 1.75 x 10" 13 s s" 1 . This is consistent with 
(but more accurate than) the incoherent analyses by Z02 
and K02, and establishes that RX J0720. 4-3125 is spinning 
down. 



4 PERIOD DERIVATIVE FROM THE 
INDIVIDUAL DATASETS 

The new data allow a more accurate determination of the 
period derivative using a conventional least squares fit to the 
individual periods than was possible in Z02. This is because 
the quality of the recent data permit a breaking of the ambi- 
guities that were present in the dataset used in Z02 caused 
by the alias patterns in the ROSAT data. We perform a 
new analysis, both in order to provide a robust estimate of 
P which in itself is now sufficiently accurate to provide sci- 
entific constraints, and in order to provide a defined region 
of parameter space for the more computationally intensive 
phase- coherent analysis which follows. 

Given the ad-hoc nature of the epochs at which ob- 
servations have been taken, it is necessary to proceed in a 



5 PERIOD DERIVATIVE FROM A 
COHERENT ANALYSIS 

The change to the model used in the MLP noted above, to- 
gether with the revised ROSAT timings (unimportant in the 
incoherent analysis above) means that the coherent period 
analysis in section 4.3 of Z02 must be revised - and indeed 
the 99 percent confidence limits in Figure |3] exclude all of 
the four identified Po, P pairs in their table 3. 



5.1 Individual Satellites 

As explained in Z02, the time taken to compute coherent 
2-dimensional MLPs for a dataset of this size is prohibitive, 



Timing analysis of RX J0720.4-3125: II 5 





1 r ■ - 







chi2_lir8a.ml.ps (2003 Jul 21 at 09:14:38) 



l.OE+08 -5.0 



since R93 start (sec) 












-l.OE+08 -5.0E+07 O.0E+00 S.OE+07 



l.OE+08 l.SE+08 2.0E+08 2.5E+08 
le since R93 start (sec) 




Figure 2. The periods as derived from the MLPs of the XMM- 
Newton, ChOl, R96d and R93 datasets (Tabled shown as a func- 
tion of time since the start of the R93 data. Error bars show la 
uncertainties. The solid line is the best linear fit with dotted lines 
showing the 99% confidence interval. The period aliases for the 
R98 and ChOO data are also shown. Only a single alias peak is 
admissable in each case (top). The best linear fit including the 
selected R98 and ChOO aliases. Also shown are the period aliases 
from the combined R96a and R96b data. Only a single alias peak 
is admissable (centre). The best fit including the selected aliases 
from the R96a and R96b data. Also shown are the period aliases 
from the combined R96c and R96d data. Again, only a single alias 
peak is admissable within the 99% confidence range (bottom). 



so that it has to be performed in stages. This can be done in 
several ways. Here we have analysed the ROSAT, Chandra 
and XMM-Newton datasets separately to start with. This 
avoids any residual difficulties in relating the different time 
systems used for the different satellites, or from any different 
assumptions in the data reduction software. 




3.39106 B. 39 1 07 



3.39108 H. 39109 8.39110 8.391 
Period (s) 



Figure 3. The 99%, 95%, 90% and 68% confidence interval con- 
tours in the [Po,P] plane for the data in Figured ('bottom'), and 
including the selected alias from the R96c and R96d combined 
data. The central contour is drawn only to identify the minimum 
in the \ 2 plane. 



Dataset 



Period 

(s) 



ToA 
(MJD) 



RASS 

R93 

R94 

R96a 

R96b 

R96a+b 

R96c 

R96d 

R96c+d 

S97 

R98 

ChOO 

XOOa pn 

XOOa MOSl 

XOOa MOS2 

XOOb PN 

XOOb MOSl 

XOOb MOS2 

ChOl 

X02 pn 

X03 pn 



8.390935±0. 
8.39120 ±0. 
8.39065 ±0. 
8.3913 ±0. 
8.3925 ±0. 
8.391100±0. 
8.392 ±0. 
8.391137±0. 
8.391095±0. 
8.39113 ±0. 
8.391109±0. 
8.391121±0. 
8.391110±0. 
8.391110±0. 
8.391113±0. 
8.391122±0. 
8.391060±0. 
8.390850±0. 
8.391115±0. 
8.391124±0. 
8.391136±0. 



000055 

00045 

00005 

0008 

0015 

000007 

025 

000063 

000005 

00012 

000013 

000020 

000014 

000047 

000034 

000050 

000170 

000190 

000011 

000002 

000015 



48176. 
49257. 
49484. 
50198. 
50210. 
50204. 
50353. 
50391. 
50372. 
50523. 
50925. 
51577. 
51677. 
51677. 
51677. 
51869. 
51869. 
51869. 
52248. 
52585. 
52761. 



0202551±0. 
2547035±0. 
0535939±0. 
6873653±0. 
5563095±0. 
6220263±0. 
9975527±0. 
3007509±0. 
8348286±0. 
7056342±0. 
6882151±0. 
0395669±0. 
4067206±0. 
4688770±0. 
4712089±0. 
9571032±0. 
9949797±0. 
9949796±0. 
6768196±0. 
9688278±0. 
9950590±0. 



0000059 
0000028 
0000092 
0000031 
0000058 
0000039 
0000064 
0000016 
0000017 
0000040 
0000036 
0000026 
0000003 
0000010 
0000007 
0000006 
0000012 
0000014 
0000008 
0000003 
0000004 



Table 2. The periods and ToAs derived from the individual 
datasets. 



In each case we have searched the [-Po,-P] plane in the 
region within the 99% confidence intervals in Figure [3] using 
the 2-dimensional MLP as described in Z02 and modified as 
described in Section 3 above. Care was taken to ensure that 
the region was searched on a grid finer than the Nyquist 
frequencies to ensure complete sampling of the \ 2 plane. 
In the case of the ROSAT data, we have included the 1993, 
1994 and 1996 data, with a total timespan of 3 years, similar 
to that for the XMM-Newton data, while the two Chandra 
observations span ~ 20 months. This resulted in MLPs with 
dimensions of 2000 x 2000, with datasets of 1.4 x 10 6 , 8 x 10 4 
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and 2 x 10 4 events in the case of XMM-Newton, Chandra and 
ROSAT respectively. 

The results are shown in Figure 0] The outer contours 
are the 99% confidence intervals for two parameters of in- 
terest, except for the XMM-Newton MLP where contours 
to include 99.9% significance are shown for illustration. The 
Chandra MLP shows diagonal alias patterns resulting from 
the availability of data at only two epochs. The ROSAT 
MLP consists of two main diagonals, broken into regions 
of better and less good fit, together with a lower likelihood 
contour parallel to these at shorter periods. In the case of 
the XMM-Newton MLP there are five regions, although at 
the 99% confidence level, only the central upper left region 
is significant. 

Overlaying these MLPs, it is clear that there is no sin- 
gle location at which all three datasets intersect. The con- 
sequence is that there is an inconsistency between the three 
determinations. In particular, the [Po,P] contours from the 
XMM-Newton MLP do not intersect with those of either the 
ROSAT or Chandra MLPs. This has persisted, despite an 
exhaustive reassessment of the analysis. 

5.2 Simulated data 

In order to investigate the viability of the coherent analysis 
(K02 questioned the applicability of approach in the study 
by Z02) and to verify more stringently the analysis process, 
we have generated simulated data with the same mean lev- 
els, noise characteristics and observation durations as the 
real data in Table 1. For those observations which were bro- 
ken into several shorter sections (such as the Chandra and 
the ROSAT data), the simulated data were generated sim- 
ilarly. An intensity variation of linearly increasing period, 
with Po and P similar to that deduced from the incoherent 
analysis was introduced into the data. The simulated dataset 
therefore closely matched the real dataset, with the differ- 
ence that the Po and P were known, and the time reference 
for all the data was known with certainty to be common. 
These data were phase folded using equation (4) and the 
known Po,P, and it was checked that all data co-phased 
correctly. 

We proceeded with the analysis as follows. The [Po,P] 
plane bounded by [8.391050 - 8.391120, - 2 x 10~ 13 ] was 
searched using the MLP, first coherently within the ChOO 
and X00a,b datasets, then in the ChOl and X02a,b datasets 
and finally in the R96 datasets. Each search resulted in a 
diagonal set of contours in A% 2 , corresponding to the alias 
patterns. We verified that the known [Po,-P] lay on one of 
the best-fit valleys. The orientation of the contours is in- 
creasingly angled with time, so that overlaying the three 
2-dimensional MLPs identified a limited number of loca- 
tions where the three datasets gave best fits. We confirmed 
this more quantitatively by a simple multiplication of the 
three MLPs. One of these loci (but not quite the minimum) 
was the input [Po , P] ■ We then carried out a coherent MLP 
analysis of all of the R96, ChOO, X00a,b, ChOl, and X02a,b 
datasets in the limited [Po, P] spaces around these loci. This 
coherent analysis generated the A^ 2 values requiring mu- 
tual co-phasing of all of the datasets, and now the minimum 
in the plane was found indeed to coincide with the input 
[Po,P]. 

This exercise verified both the MLP analysis routines 
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Figure 4. Contours of A\ 2 in the [Po, P] plane for the R93, R94 
and R96 dataset (top) for the ChOO and ChOl dataset (centre) and 
the X00, X02 and X03 dataset (bottom). In the XMM-Newton 
data, were the 99% confidence interval contours (rather than the 
99.9% in this case) only to be shown, then only the central region 
at the top left is significant. 
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in as realistic fashion as possible, and the feasibility of iden- 
tifying the correct period amongst the alias patterns, given 
the nature of the particular datasets at hand. 

We repeated this order of analysis for the real datasets. 
Overlaying the MLPs identifies three loci where there is 
overlap in the patterns, around [8.3910771, 1.796 x 10" 13 ], 
[8.3910954, 9.50 x 10~ 14 ] and [8.3910865, 1.3825 x 10" 13 ]. A 
detailed search using the entire dataset identifies the region 
around the first solution as the only region where a coher- 
ent phase can be maintained over the entire dataset. The 
best fitting [Pq,P] = [8.39107712, 1.7956 x 10~ 13 ]. However, 
when folded, this [Po,P] pair is inconsistent with the X03 
dataset. 

5.3 Possible origins for the discrepancy 

At this stage we are not able to arrive at a consistent co- 
herent solution for a [Po, P] pair from our extensive dataset. 
However: firstly, we have shown from the simulations above 
that the dataset in principle contains sufficient information 
to identify the correct solution, even if we do not include 
data from the X03 observations, or from S97, R93 or RASS 
datasets. Therefore, if a solution with P constant does ex- 
ist, the technique based on the coherent analysis allows us to 
identify it. Secondly, we have also verified that the software 
used to calculate the MLP returns the input parameters for 
the simulation correctly: we believe that any coding errors 
are therefore unlikely. Thirdly, in our analysis in Section f5. II 
we have taken care not to mix the data from different satel- 
lites, so we have avoided any residual difficulties in relating 
the different time systems used for the different satellites. 

Consequently we conclude that either the event tim- 
ings used in our data are incorrectly assigned with respect 
to observations at other epochs with the same satellite, or 
that the model we have used for a constant period change 
[equation (5)] is inappropriate. We discuss these in turn. 

We have cross-checked the event timings for our data 
where they overlap with K02. K02 have revised their table 1 
post-publication, in the version held on the Astro-ph server 
which corrects the time system reference and other errors 
used to calculate their ToAs in their published version. We 
have carried out a detailed cross-check between the ToAs in 
their revised table and our datasets in Table 1. Our ToAs 
are given in Table [5] For those datasets in common (R93, 
R96d, S97, ChOO and ChOl), the ToAs we calculate match 
those derived by K02 closely (to within 2% of the period), 
and we derive similar uncertainties. This indicates that these 
independently-derived datasets have been reduced consis- 
tently, and the likelihood of systematic errors remaining is 
now small. However, it is impossible to eliminate the possi- 
bility of error entirely since similar misconceptions could be 
resident in the software, or use of the software, for reducing 
the data - such as was the cause of the error in relating 
ROSAT and Chandra data in Z02 and K02. 

In the case of the XMM-Newton timings, no such cross- 
check is currently possible. It is also of some concern that 
the best fit [Po, P] pair from the XMM-Newton observations 
alone (top centre contours in Figure^ are located outside of 
the 90% confidence intervals calculated from the incoherent 
analysis (Figure |3J. We have therefore investigated the tim- 
ing accuracy of XMM-Newton in further detail. After the 
XMM-Newton SAS barycen task was corrected for version 



13.1, an absolute time shift in Crab pulsar data of ~ 1.2 
msec remained with respect to the radio data. This resulted 
from an incorrect sign used to compensate for internal elec- 
tronic delays. Becker (private communication) reports that 
there is still a time difference of between 600 and 300 fisec in 
the Crab data (apparently decreasing with time) so minor 
timing discrepancies remain. 

These discrepancies appear to be on such a scale that 
they would not affect our XMM-Newton analysis. How- 
ever, there is some indication that there may be larger 
scale discrepancies inherent in some XMM-Newton datasets. 
These have arisen in some eclipsing binary data, for example 
OY Car (Wheatley & West 2003), and in EP Dra (Bridge 
2004). In the former case, the eclipse occurs 55 sec early, 
while in the latter it occurs 69 sec late. These residuals 
are significantly larger than permitted by the uncertainties 
in the ephemerides for the two stars and would invalidate 
the coherent analysis for our XMM-Newton datasets of RX 
J0720.4-3125. Schwope et al. (2004) also find some timing 
discrepancies for DP Leo, but report on the other hand that 
the eclipses in XMM-Newton observations of HU Aqr are 
consistent with the orbital ephemeris. The situation remains 
unclear. 

Regarding the appropriateness of a constant P model 
[equation (5)], we note three possibilities: glitching may have 
occurred, or there may be a deviation from a constant spin 
down, or a modulation in arrival times due to the star be- 
ing in a binary system. Glitching behaviour has been ob- 
served in neutron stars with spin periods as long as that of 
RX J0720. 4-3125, as for example in the anomalous X-ray 
pulsars (AXP) 1E2259. 1+586 and 1RXS 1708-4009. Long 
term monitoring of several AXPs with the Rossi X-ray Tim- 
ing Explorer (Kaspi et al. 2001, Gotthelf et al. 2002, Gavriil 
& Kaspi 2002) revealed the diversity in the behavior of the 
single objects, ranging from high stability (in 1E2259. 1+586 
in quiescence, and 4U 0142+61 for which a linear fit with 
constant v phase-connects data collected over more than 4 
years) to instabilities so severe that phase-coherent timing 
is not possible (as in IE 1048.1-5937). In the case of AXPs, 
timing stability decreases with increasing v and the fre- 
quency derivative of RX J0720. 4-3125 is one order of mag- 
nitude lower than that of the most stable AXPs. However, 
since the two sources belong to different classes of pulsars, 
it is not really obvious how much the extrapolation of this 
trend can be trusted. 

So far, regrettably the sampling of the current dataset 
is not likely to be sufficient to distinguish whether either of 
these possibilities is the explanation for the inconsistencies 
we obtain in our MLPs. 

5.4 Combined analysis 

As discussed in § 15.11 we cannot find a consistent single 
[Po , P] pair from the analysis of the datasets from different 
satellites. However, there is an overlap between the ROSAT 
and Chandra MLPs, therefore it is possible to compute a 
coherent MLP by restricting the analysis only to these data. 
The agreement in ToAs computed by K02 and ourselves here 
suggests that errors in the event timings are unlikely. Owing 
to the sparseness of the sampling, six aliases remain in the 
coherent analysis. These all have 1.1 x 10~ 13 < P < 1.65 x 
10" 13 and are shown in Figure EJa). 
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Figure 5. The final 99% confidence contours for the coherent 
analysis of the combined R93, R94, R96, ChOO and ChOl datasets 
(top) and those for the coherent analysis of the XMM-Newton 
datasets in more detail (bottom). In the bottom panel, the ap- 
parent aliases are the result of the contouring software. 



The XMM-Newton dataset is the largest and provides 
highest S/N ratio MLPs. If we take it at face value, a region 
with P ~ 1.95 x 10" 13 is identified in the MLR This is 
shown in Figure [^Jb). 



5.5 The incoherent analysis revisited 

In the light of the above uncertainties, it is necessary to 
examine whether the incoherent analysis presented in Sec- 
tion 4 has been compromised. We consider this unlikely, be- 
cause the precision derived from a coherent analysis of the 
time spanned by our data is nearly two orders of magnitude 
greater than that which can be derived from the incoherent 
analysis. The incoherent analysis is therefore unlikely to be 
affected by the uncertainties at the accuracy required for it. 
The incoherent analysis is also immune from many sources 
of error which might come to light only when data from 
different epochs are combined, such as might arise for light 



travel time corrections, or the transformation between time 
systems. 



6 DISCUSSION 

6.1 Interpretation of the spin down rate 

We have presented a revision of the measure of the RX 
J0720. 4-3125 spin down rate published last year by Z02. 
By using additional Chandra and XMM-Newton datasets, 
we have been able to determine a positive spin down rate of 
P = (1.4 ± 0.6) x 10~ 13 s/s at 99% confidence level, as de- 
rived from the incoherent analysis. More refined values have 
been derived from a coherent analysis, separating ROSAT 
and Chandra datasets from the XMM-Newton ones. How- 
ever, a self-consistent phase-coherent solution based on the 
entire set of data from the three different satellites cannot 
be found. Due to these unresolved issues, in order to discuss 
the physical implications of our measure we concentrate on 
the result obtained from the incoherent analysis. Because 
of the larger numbers of datasets available since Z02 and 
K02, this measure of the spin-down rate is now sufficiently 
constraining for our needs. 

Despite the fact that the value of P reported here is 
higher than in Z02, the two solutions are sufficiently simi- 
lar that most of the scientific conclusions published by Z02 
remain unchanged. 

First of all, our measure rules out the possibility that 
the source is accreting from the interstellar medium, since 
this requires P < 10" 15 s/s (see Z02, K02). 

Furthermore, makes it unlikely that RX J0720.4-3125 
is spinning down under propeller torques exerted by the in- 
terstellar medium. In Z02 our argument against this was 
based on the combination of two constraints between mag- 
netic field, density of the medium and star velocity. The first 
was obtained from the expression of the propeller spin down, 
and the second by imposing that the star has entered the 
propeller phase. The latter gives B12 < 25y / nn M) 3//2 , where n 
is the external density in cm" 3 , B12 = B/(10 12 G) and vio 
is the star's velocity normalized to 10 km/s. For the former 
condition we used a propeller model as in Colpi et al. (1998) 
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that, for P = (1.4 ± 0.6) x 10~ 13 s/s, gives B 12 « (36 - 
139)n _9//8 «^Q . Therefore, the two conditions taken to- 
gether limit the star's velocity to Vio < n 1//3 (n ~ 1 for the 
interstellar medium). This modest velocity was regarded as 
unlikely at the time of the Z02 paper, and it is now ruled 
out by the recent measure of a proper motion of 97 ± 12 mas 
per year (Motch, Zavlin & Haberl, 2003) , which corresponds 
to a transverse velocity Vt ~ 50 x (D/100pc) km/s. Since 
Vio is likely to be greater than unity, the propeller mecha- 
nism dominates only if the star is presently going through a 
high density medium, with n > 100. The possibility that RX 
J0720. 4-3125 is passing through one of these high density 
clouds (which are suggested to be present on small scales 
ranging from 10 to 10 6 AU, see Lauroesch & Meyer 2000) 
has been considered by Motch et al. (2003). However, their 
conclusion is that this scenario is problematic since it may 
imply changes in the X-ray flux and blackbody tempera- 
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Figure 6. The location of 9 isolated neutron stars in the tem- 
perature/age diagram for two different models of proton super- 
fluidity (lp and 2p) and different masses (figure adapted from 
Yakovlev et al. 2003). In both panels, the diamond and the trian- 
gle mark the position of RX J0720.4-3f 25 as obtained assuming 
T°° = 0.4 X 10 B K and T°° ~ 0.9 X fO 6 K, respectively (see text 
for details) 

ture (not observed so far) and it cannot account for a large 
P > 10" 14 s/s. 

We are then left with two possible explanations for the 
cause of the observed spin down rate: propeller torque ex- 
erted by a fossil disk or spin down due to emission of mag- 
netic dipole radiation. Although the first cannot be com- 
pletely ruled out by our measure, and may better explain a 
high level of timing noise which is typical of accretors (Kaspi 
et al. 2001), we regard it as extremely unlikely. In fact, it 
requires two separate mechanisms to explain the observed 
X-ray luminosity and spin down rate (see Z02). Also, as dis- 
cussed by Perna et al. (2000) (see also Kaplan et al. 2003), 
disk models fail in modeling the spectrum of the optical 
counterpart. 

By interpreting the spin down rate in terms of mag- 
netic radiation losses and assuming a dipolar field we can 
constrain the spin down age, t s d = P/ (2P) ~ (0.6 — 
1.5) x 10 6 yr, and the present value of the magnetic field, 
B = (2.8 - 4.2) x 10 13 G. 2 The spin down age is consistent 
with the flight time of 10 6 yrs from the two birth places 
proposed by Motch et al. (2003), the Sco OB2 complex and 
the Vela OB2 + Trumplcr 10 association. 

6.2 Location in the temperature/age diagram and 
estimate of the mass 

The age derived from our timing analysis can be then com- 
pared with the star's age as inferred by the cooling curves, 

2 As in Z02, we reiterate that this must be regarded as an upper 
limit. A twisted magnetosphere may lead to a reduction up to 
an order of magnitude in the inferred polar value of the magnetic 
field. The estimate of the spin down age is unaffected (Thompson, 
Lyutikov & Kulkarni, 2001). 



which is based on the spectral measure of the surface tem- 
perature. 

Different theoretical cooling models have been pre- 
sented in the literature. Roughly speaking, they predict a 
two-fold behavior of the cooling curves: slow cooling for low- 
mass neutron stars and fast cooling for stars with higher 
masses. The transition between the two regimes can be more 
or less sharp, depending on the assumed superfluid proper- 
ties of the neutron star interior. It has been realized (see 
Yakovlev et al., 2003 for a recent discussion) that simple 
models which do not account for proton and neutron super- 
fluidity fail in explaining the surface temperatures observed 
in many sources, unless there is a fine tuning in the masses 
of objects such as Vela, Geminga and RX J 1856-3754, and 
all of them have nearly the critical mass that bounds the 
(in this case, sharp) transition between slow and fast cool- 
ing regimes. Models with proton superfluidity included do 
not require this unlikely assumption and predict the exis- 
tence of a region, intermediate between the two regimes and 
relatively large in the parameter space, which is populated 
by medium mass neutron stars (roughly between 1.4 and 
1.65 Mq). Although the full picture only holds if, at the 
same time, neutron superfluidity is rather weak, so far the 
latter models provide a better comparison with observations. 
It is interesting that many neutron stars (as IE 1207-52, 
Vela, RX J1856-3754, PSR 0656+14) have a surface tem- 
perature which falls in such a transition region (see Figure|nj 
re-adapted from Yakovlev et al., 2003). 

If the origin of the X-ray flux from RX J0720.4-3125 
is the thermal emission from the cooling neutron star, in 
first approximation the surface temperature can be taken 
equal to the temperature of the blackbody fit, i.e. T°° ~ 
0.9x 10 6 K (Haberl et al., 1997; Paerels et al., 2001; Haberl et 
al., 2004). This is likely to be an upper limit, since most de- 
tailed atmospheric models or fits which account for a broad- 
band deviation from a blackbody tend to give a smaller tem- 
perature (see e.g. Pavlov et al., 2002; Burwitz et al., 2003). 
As we can see from Figure HJ this gives a maximum cooling 
age of ~ 0.3 x 10 6 yr (for M ~ 1.35 Mq) which, given the 
numerous uncertainties, may be considered as marginally 
compatible with the spin-down age. Alternatively, Kaplan 
et al. (2003) propose a multicomponent fit to optical/X-ray 
data of RX J0720. 4-3125 concluding that the X-ray emis- 
sion at ~ 0.9 x 10 6 K originates from hot polar caps, while 
the whole surface is cooler with a blackbody temperature of 
T°° ~ 0.4 x 10 6 K. If this is the case, the spin down age 
is compatible with RX J0720. 4-3125 being a medium mass 
neutron star with M ~ 1.5 — 1.6 Mq, depending on the 
kind of proton superfluidity assumed in the model (lp and 
2p respectively). 

It is worth noting that given the significant spread in 
the cooling curves in Figure [5] resulting from a strong mass 
dependence, an age as low as 10 4 yr (which may be possible 
if the magnetic field decays rapidly) is also consistent with 
the observed surface temperature. We discuss these issues in 
the next section. 

6.3 Implications for the past history of RX 
J0720. 4-3125 

The considerations presented so far are based on the spin 
down age t a d, which is representative of the true age of the 
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source only in the case in which the magnetic field has re- 
mained almost constant during the star evolution. The same 
condition applies for the validity of the cooling curves men- 
tioned above, which do not include the extra input of energy 
released in the neutron star in case of field decay. Taken 
face value, the timing parameters of RX J0720. 4-3125 are 
compatible either with those of a radio-quiet cooling pul- 
sar with high (but not extreme) magnetic field or an old 
magnetar, born with B > 10 14 G and still kept warm by the 
decay of its superstrong field. As in Z02, in order to discrimi- 
nate between these possibilities, we studied the evolutionary 
tracks on the B, P diagram. Three different mechanisms are 
typically proposed for inducing field-decay: ambipolar dif- 
fusion in the solenoidal or irrotational mode and Hall cas- 
cade (Goldreich & Reisenegger 1992, Heyl & Kulkarni 1998, 
Colpi, Geppert & Page 2000). In reality all the three pro- 
cesses co-exist with different time scales, and each of them 
may dominate the evolution depending on the istantaneous 
values of B, L. In absence of more detailed computations, 
Heyl & Kulkarni (1998) and Colpi, Geppert & Page (2000) 
tentatively isolated the three processes and proposed some 
simple, phenomenological rules to mimic the evolution in 
the three regimes. As in Z02, we used as first approxima- 
tions these decay rules and re-computed the source age, Td, 
and the value of the magnetic field at the birth of the neu- 
tron star, Bo corresponding to our revised measure of P. By 
measuring Bo in 10 13 G, P in seconds, P in s/yr and the 
age in 10 6 yr, the resulting expressions are 
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where b ~ 3, Po is the period at the star birth, and the 
parameters a, a discriminate between the three decay laws 
(Colpi, Geppert & Page, 2000). The results are shown in 
Table In all cases the star is assumed to be born with a 
period of 1 ms: results are not strongly dependent on this 
exact value, provided it is less than the present period 3 . 

As we can see, the only scenario compatible with a su- 
perstrong field at the star's birth (Bo ~ 10 G) is that 
one involving a very fast decay, i.e. the Hall cascade. In this 
case the predicted star's age is ~ 4 x 10 4 yr. Although such a 
young age is only marginally compatible with the absence of 
a remnant, it is not incompatible with that inferred from the 
cooling curves (see Figure |SJ provided the mass is slightly 
larger. A slightly larger mass still is required if we allow for 
the extra-heating due to B-decay from Hall cascade (Gep- 
pert & Colpi, private communication). 

On the other hand, both mechanisms involving ambipo- 
lar diffusion predict a magnetic field which is relatively con- 
stant over the source lifetime and close to the present value. 
The corresponding star's age is between ~ 0.4 x 10 6 yr and 
« 1.4 x 10 6 yr, and is close to r s d. 



3 K02 discuss the possibility that RX J0720. 4-3125 is an example 
of the "injection" hypothesys (Vivekanand & Narayan, 1981) and 



Solution 


B-Decay Mechanism 


Bo 


age 






10 13 G 


(years) 


1 


Hall Cascade 


120.2 


3.5 x 10 4 


1 


Amb. diff . , irrotational mode 


3.0 


1.4 x 10 6 


1 


Amb. diff., solenoidal mode 


5.3 


0.8 x 10 6 


2 


Hall Cascade 


121.6 


2.3 x 10 4 


2 


Amb. diff., irrotational mode 


4.4 


0.6 x 10 6 


2 


Amb. diff., solenoidal mode 


7.0 


0.4 x 10 6 



is born with a long period, Po ~ 
laws predict an age of order 10 4 



3.3 s. In this case, all three decay 

yr- 



Table 3. Predicted source age and primordial field for three dif- 
ferent mechanisms of decay, simulated as in Colpi ct al. (2000). 
Here P = 8.391 s, and we used different values of P: P = 
0.8 X 10~ 14 s/s (solution '1') and P = 2 X 10~ 13 s/s (solution 
'2'). These correspond to the boundaries of the 99% confidence 
range derived from our coherent analysis. In all cases, the source 
is assumed to be born with P = 1 ms. 



7 CONCLUSIONS 

We have eliminated the discrepancies between the calcu- 
lated ToAs in the K02 paper and those derived from the 
same data analysed by Z02. This invalidated our derived 
[Po, P] solutions in Z02 (K02 have also updated their ToAs). 
In order to correct these solutions, we have reanalysed the 
datasets and taken the opportunity to incorporate signifi- 
cant new datasets from both Chandra and XMM-Newton. 
At the same time, we have revised the model used by Z02 to 
be fitted to the data to retain the correct phase reference. 

Taken individually, the combination of old and new 
datasets allow a significantly more accurate [Po,P] solution 
to be derived than that in K02 or Z02. This establishes that 
RX J 0720. 4-3125 is spinning down. 

We have also attempted a phase-coherent analysis of 
the combined datasets, as in Z02. However, we have been 
unable to find a consistent solution, even when we analyse 
the data from different satellites separately. In this case the 
XMM-Newton data are not consistent with the other data. 
We have explored the possible reasons for this: they may be 
intrinsic to RX J0720. 4-3125 (in which case the constant P 
model is not appropriate) , or they may result from continu- 
ing timing errors in the data we have not succeeded in iden- 
tifying. We present the solutions from the combined ROSAT 
and Chandra data and the XMM-Newton data separately. 

The individual (incoherent) [Po,P] solution is suffi- 
ciently accurate to allows us to place strong constraints 
on the mechanism causing the spin-down. The large P ~ 
10~ 13 s/s rules out accretion and propeller origins and points 
toward an interpretation in terms of magneto-dipolar break- 
ing. This in turn constrains the age of RX J0720. 4-3125 and 
its magnetic field strength. 

For a magnetic field as strong as B = (2.8—4.2) x 10 13 G, 
the electron cyclotron line is unaccessible (it falls at about 
350 keV). However, a relatively broad proton cyclotron ab- 
sorption feature is predicted at ~ 0.2 — 0.3 keV (Zane et 
al., 2001). The presence of a cyclotron absorption line at the 
low energy of the sensitive energy band of the XMM-Newton 
EPIC instruments has been indeed recently reported (see 
Haberl et al., 2003a, 2004). The feature may explain the ob- 
served deviations from a Planckian shape at these energies, 
and possibly also the anticorrelation of the modulation of 



Timing analysts of RX J0720. 4-3125: II 11 



hardness ratio and total X-ray intensity detected by XMM- 
Newton (Cropper et al., 2001). 

Our study allows a possible evolutionary link between 
dim INS and AXPs (which relies on the interpretation of the 
former as warm-out magnetars) but only if the INS have an 
age of ~ 10 4 yr. A more conservative interpretation is that 
RX J0720. 4-3125 was born with a strong, but not super- 
strong, field [B = (2.8 - 4.2) x 10 13 G] which is compatible 
with those of the canonical radio-pulsars. Similar conclu- 
sions have been argued also by K02, however their measure 
of P has a large error and does not permit the discrimination 
between this and the separate possibilities that RX J0720.4- 
3125 has a more conventional field (B ~ 10 12 G) and lower 
spin-down rate (P ~ 10 _lj s/s) or is even spinning-up. Com- 
paring with the radio-active pulsar population, objects with 
B > 10 13 G are rare, but their evidence is rapidly grow- 
ing (Camilo et al., 2000). The parameters of RX J0720.4- 
3125 we have derived are not too dissimilar from those of 
PSR J1814-1744, which has P ~ 4 s and P ~ 7.4x 10~ 13 s/s. 

One of the most striking mysteries about ROSAT iso- 
lated neutron stars, AXPs, soft-7 repeaters and objects 
like Geminga is why they do not exhibit radio emission. 
As suggested by Motch et al. (2003), who derived very 
similar conclusions basing on the proper motion measure- 
ment, it could thus be that RX J0720. 4-3125 and several 
others apparently radio-quiet neutron stars are just radio 
pulsars with off-axis beam that does not cross the Earth. 
The radio beam narrows with increasing periods (Briggs 
1990), making this explanation even more plausible. On the 
other hand, evidence for a population of genuinely radio- 
silent young neutron stars arises from population synthe- 
sis (Neuhauser & Triimper, 1999; Gotthelf & Vasisht, 2000; 
Popov et al., 2000a, b and references therein). The group 
of ROSAT INSs is very homogeneous and they are all rela- 
tively close-by (within 300-400 pc) , making unplausible that 
all their beams are misaligned. At present, these two possi- 
bilities cannot be distinguished. Detailed timing analysis of 
other ROSAT INSs is crucial to reach a more comprehensive 
understanding of the entire population. 
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